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Abstract 

We study the prospects of probing large extra dimension model at the LHC 
through neutral triple gauge boson production processes. In theories with 
extra dimensions these processes result from the exchange of a tower of mas- 
sive graviton modes between the SM particles. We consider 777, 77.Z', 7ZZ 
and ZZZ production processes, and present our results for various kinematic 
distributions at the LHC for y/S = lA TeV. 
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1 Introduction 



The results that would emerge from the Large Hadron Collider (LHC) are expected 
to reveal patterns under which the standard model (SM) is likely to be modified 
at the TeV scale. At these high energies, LHC is expected to unveil the mystery 
of large hierarchy between the electroweak and the Planck scales and also produce 
interesting new physics signals. Two well-studied scenarios of physics beyond the 
SM at the TeV scale are extra dimension theories and supersymmetry, which will be 
tested at the LHC among many others. The large extra dimension model, proposed 
by Arkani-Hamed, Dimopoulos and Dvali (ADD) [1], circumvents the hierarchy 
problem by exploiting the geometry of extra spatial dimensions to find a mechanism 
to lower the Planck scale. The tower of Kaluza-Klein (KK) modes as a result of 
gravity propagating the full 4 + dimensions would mediate scattering processes and 
hence lead to deviations from the SM predictions. At colliders, exchange of virtual 
KK modes or emission of real KK modes give rise to interesting phenomenological 
signals at the TeV scale [21 [3] . Virtual effects of KK modes lead to the enhancement 
of the cross section of pair production in processes viz. di-lepton, di-gauge boson 
(77, ZZ, W~^W~), dijet [H El El El E] in some kinematic regions. The real emission 
of KK modes leads to large missing signals viz. mono jet, mono photon, mono 
Z boson and mono boson [9]. 

The di-gauge boson final states have been extensively studied in the context of 
extra dimension models. The triple gauge boson final state is an interesting new 
physics signal in some of the beyond SM scenarios [TU] . In this paper we consider the 
neutral triple gauge boson production at the LHC and study how the ADD model 
would alter the SM expectation. In the SM, the triple gauge boson final state is 
an important signal as it depends on the 3-point and 4-point couplings among the 
gauge bosons which is a test of the electroweak theory. This process in the SM has 
been studied to leading order (LO) [HI [T2] and its extension to the next to leading 
order (NLO) was on the Les Houches wish list [10] and has been finally achieved 
recently [131 [HI [151 [IS]- The triple gauge boson production processes in the SM 
are the precise predictions of the electroweak gauge theory and gauge self-couplings. 
They are also potential background to many new physics models like supersymmetry 
and Technicolor. For example, ^'77 in SM is a background to signals with di- 
photons and missing transverse energy in gauge mediated super symmetric theories 
[T7] and 777 production in SM is a background to one photon plus techni-pion [18] . 
Processes with three gauge bosons can also come from the ADD model as gravitons 
couple directly to gauge bosons of SM. While mono-jet or di-lepton production 
is more sensitive to parameters of models with extra-dimensions compared to the 
triple gauge boson production, all these processes involve same universal coupling of 
gravity with the SM particles, and hence can provide equally important information 
about the model. Moreover, in discriminating physics beyond the SM namely SUSY 
or technicolour models using triple gauge boson production, one can not ignore the 
potential contributions resulting from models with extra dimensions. 
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In this analysis we consider the process PP — )■ VVV X, where we restrict to the 
neutral gauge bosons \^ = 7, Z and X is some hadronic final state. The following 
four final states are the subject of this analysis: (i) 777 (ii) 77Z (iii) 7ZZ and (iv) 
ZZZ. The case where V = will be part of a different paper [19]. 

This paper is organised as follows: in section 2 we present the analytical calcula- 
tion of the above mentioned processes with a brief introduction to the ADD model, 
section 3 is arranged to present the numerical results of our studies and finally we 
summarise the results in the last section. 



2 Neutral triple gauge boson production 

The TeV scale colliders like the LHC are suitable for studying the quantum gravity 
effects in the ADD model. This model is defined in such a way that the SM particles 
are localised on a 3-brane with negligible tension and only gravity can propagate 
in the full 4 + d dimensions as it describes the geometry of the full space time. 
The extra spatial dimensions d can be compactified on a d-dimensional torus with a 
common radius R/ (27r) of macroscopic size. As a consequence, the effective Planck 
scale (Mp) in 4 space-time dimensions is related to the fundamental Planck scale 
(Ms) in 4 + d dimensions as follows, 

M| = Cd Ml+^ , (1) 

where Cd = 2(47r)~^/r((i/2). Although, the SM particles are localised on a 3- 
brane, they can sense the extra dimensions by their interaction with gravity via the 
massive KK towers of the graviton. The corresponding interaction Lagrangian is 
the following, 

00 

^ = -^E^'^(^)^J(^) ' (2) 
^=0 



where k = y/lGn/Mp, T^"" is the energy-momentum tensor of localised SM fields 
and hj^^ denotes the massive KK graviton labeled by a d-dimensional vector it with 
all positive components, hj^^^ contains one spin-2 state, (n — 1) spin-1 states and 
n(n — l)/2 spin-0 states, having mass. 

The zero mode of the KK tower corresponds to the massless graviton in the 4 space- 
time dimensions. The effective graviton propagator, after summing over all KK 
states, can be expressed as. 



m^. + ie 

rl 
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(^j [tt + 2^/(A/y^)] , (4) 
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where s is the invariant mass of the boson pair resulting from the decay of the KK 
mode and the function ^/s) is described in pjj, which depends on the ultra violet 
cut-off A. Although the interaction of KK modes with the SM particles is suppressed 
by the coupling n (Eq. ([2])), the cumulative effect of summing over large number of 
accessible KK modes (Eq. (jl])) compensates the suppression, making the effective 
coupling significant enough to have observable effects. It is usual practice to set the 
UV cutoff A = Ms and simplify the summation of virtual KK modes [21 [3] to do 
the phenomenology. We follow the approach of [3] for this analysis which retains 
the details of the number of extra dimensions. 

The neutral gauge boson final state at the hadron collider PP — t- VVV X at 
LO comes from the sub process 

q{pi) + (i{p2) v{p3) + v{pi) + v{p5) , (5) 

where = 7, Z and X is any final state hadron. The SM diagram for the above 
process is shown in Fig. [1] with all possible permutations of final states. For the final 
state with at least two ZZs, the Higgs could contribute by coupling to the quarks, 
but this is negligible in the vanishing quark mass limit. In the case of ZZZ final 
state, there are additional Higgs strahlung diagrams but their contribution is also 
quite small and vanishes in the large Higgs mass limit. Hence, we have not included 
the processes with the Higgs boson. In the ADD model, the KK modes couple to 

, [WWWW 

WWWW 'If) 
9W , vvwww 

Figure 1: Typical Feynman diagram for triple gauge boson production in SM. 

V bosons, quarks, anti-quarks as well as to quark-antiquark-V boson vertex. The 
four categories of Feynman diagrams that give a VVV final state in ADD model 
are shown in Fig. |2l We have used unitary gauge — )■ 00) for Z-boson and the 
Feynman gauge = 1) for the photon. 

In the SM, the LO process for the production of 777 at hadron colliders results 
from the annihilation of a quark and an anti-quark. In the ADD model, the produc- 
tion mechanism is again from the same initial states but one of the photon comes 
from either quark or an antiquark and other two photons come from the decay of KK 
graviton. The typical Feynman diagrams that contribute in the SM and in the ADD 
model are shown in Fig. [H [2j The Feynman rules for the processes with KK gravi- 
ton can be found in [21 13] . All the expressions for the matrix element squared with 
proper spin, color sums and averages are obtained using a symbolic program based 
on FORM [20]. The KK graviton propagator in the ADD model is proportional to 
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Figure 2: Typical Feynman diagrams for triple gauge boson production in ADD 
model. Dashed line represents the KK mode of the graviton and the other particle 
lines are same as in Fig. 1. 



V 
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(6) 



n 



■ij - mh. + is 



where the invariants are Sij = {pi +PjY, and V can be evaluated using Eq. (j4]). The 
numerator of the spin-2 propagator [3j of the KK graviton is given by 



where 77^;^ = g^^ — k^ky/m^. Here k is the momentum flowing through the prop- 
agator. Terms proportional to negative powers of mass of KK mode in rj^y do not 
contribute as they are proportional to k^k^. This provides a useful check on our 
calculation. The matrix elements have been checked for gauge invariance. We per- 
formed similar computation for evaluating the parton level subprocesses for 77Z, 
'-yZZ and ZZZ productions. In the following we list few of the important observa- 
tions. 

For the 77^' production, in the limit mz — )■ 0, we reproduce the matrix ele- 
ments for 777 process with the changes: {Cy + C^)/4 — )■ Qj, — )■ e, where Cy, 
Ca are the vector and axial vector couplings of the weak gauge boson respectively, 
Tz = e/{sm9u, cos9w) and Qf is the electric charge of the quark flavors. In the case 
of 'fZZ production, we find that the parton level subprocesses in SM and ADD 
model are similar to those of the 77^' production with the changes 7 ^ Z. The 
squared matrix elements for '-/ZZ production that come from ADD model alone are 
not related to those of 777 production. The reason is that some of the terms pro- 
portional to m| that appear in the graviton-ZZ vertex cancel all the inverse power 
of m| present in the Z boson polarisation sum, giving contributions that have no 
analogous ones in the 777 process. However, the expression for SM squared matrix 
elements of 7ZZ is related to that of 777 process in the SM if we take mz — > 0, 
{Cy + QCyC\ + C^)/16 — )■ Qj and — > e. For ZZZ production, squared ma- 
trix elements involving ADD vertices do not have any relation with those of 777 



{k) 
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production for the same reason as production case. The SM squared matrix ele- 
ments of this process are related to those for the 777 process in SM with the following 
replacement @ in the limit -> 0, (C^ + IbC^C^ + 15C^C\ + C^)/64 — > Qj, 
Tz — > e. The expressions for the squared matrix elements of the processes discussed 
above are too large to be presented in this paper. 

3 Numerical results 

In this section, we present different kinematical distributions for the production of 
neutral triple gauge bosons. The predictions are for the LHC at center of mass energy 
\/S = 14 TeV. We have used CTEQ6L parton densities [21]. For the strong coupling 
constant that appears in CTEQ6L, we use Aqcd = 0.226 GeV and Uf = 5 flavors. 
We set the factorisation scale fip = Pt the transverse momentum distribution 
of V and = Q for the invariant mass (Q) distribution of the di-boson pair. In 
addition we apply the following cuts on and the rapidity y^: 

P^'^ > 25 GeV and y^'^ < 2.7 . (8) 

We also ensure that in general the invariant mass of the di-boson {i.e. any two 
identical bosons among the 3V) is less than M5. We use rriz = 91.1876 GeV and 
sin^6'^ = 0.2312. The fine structure constant is taken as a = 1/128. 

Recently, both CMS [22j and ATLAS [23] reported, searches for signatures of 
extra dimensions in the di-photon mass spectrum at the LHC for 7 TeV p p collisions. 
The 95 % lower bound on Ms vary between 2.27 - 3.53 TeV depending on the 
number of extra dimensions c? = 7 — 3 for ATLAS and Ms vary between 2.3 - 3.8 
TeV depending on the number of extra dimensions = 7 — 2 for CMS, both using a 
fixed K factor of about 1.6 [6j. Hence we have used the phenomenologically viable 
ADD model parameters Ms = 3.5 TeV and (i = 3 for present study. 

For the processes involving more than one photon, it is important to isolate 
photons from each other i.e. they need to be well separated in phase space so that 
they can be identified as separate objects in the detector. To do this we consider 
a cone of radius R = ■>/ (Ai/)^ + (Af/))^ in the rapidity-azimuthal angle plane {y, (j)) 
and ensure that the minimum separation between any two photons is taken to be 
Rryry = 0.4. lu thc foUowlug, we describe our findings for the various triple gauge 
boson production processes. 

1. 777 production: In this case, the three photons in the final states are classified 
in such a way that P^^ > P^^ > P^^. We have compared our predictions for P^^ 
distribution in the SM against those given in [16] and found a very good agreement 
confirming the correct implementation of our analytical results in our numerical 
code. In the left panel of Fig. |3l we present the transverse momentum distributions 

^We find that the most general result for the replacement of n number of Z-boson with photon 
in the SM squared matrix element is icl+clr+2n{n~iy,clcl{cl+clr-^) — ^ g2n^ ^^^^^ ^^^^^ 
for all the above three processes for n = 1, 2, 3. 



5 



T)-anverse Momentum Distribution of y(1) Scale Variation of Rapidity Distribution of y(1) 

200 400 600 800 1,000 1,200 -2 -1 1 2 




Figure 3: Transverse momentum distribution of 71 [i.e, with maximum Pt] (left 
panel) at the LHC for Ms = 3.5 TeV and d = 3. Rapidity distribution of 71 (right 
panel) at the LHC for ADD model parameters Ms = 3.5 TeV and d = 3 in the 
region where P}^ G (750, 1250) GeV. Dependence of rapidity distribution of 71 on 
the factorisation scale in the range fip = 0.2P^^ and fip = 2P^\ 

of 7i in SM as well as in SM+ADD. We have chosen Ms = 3.5 TeV and d = 3 as 
representative parameters of the ADD model. In high P^^ region, the distribution for 
SM+ADD is fully controlled by processes coming from ADD model and is enhanced 
due to the dominant contributions of the KK modes. In the right panel of Fig. [31 
rapidity distributions of the most energetic photon 71 are shown for 750 < P}^ < 
1250 GeV in SM and SM+ADD. It is seen that the SM contribution is extremely 
small in this range. 

In order to estimate the factorisation scale fip dependence present in our LO 
results, in the right panel of Fig. [3] we have plotted rapidity distributions for three 
different choices of fip i-e-, fJ^F = (0.2, 1,2)P^\ In the central rapidity region the 
variation of the rapidity distribution with respect to the factorisation scale is the 
largest. With respect to the central choice of fip = , the variation is about 
23.6 % and 8.2 % for the choice of fip = 0.2 P^^ and fip = 2 P^^ respectively. 

The Pt distribution of 72 is found to be similar to that of 71, but is different for 
73 (the least energetic of the three photons) as shown in Fig. H] (left panel). Similarly 
its rapidity distribution Fig. H] (right panel) is also different from the most energetic 
photon. 

2. 77^' production: Here, the invariant mass distribution of the photon pair 
is a useful observable because in the ADD model the photon pair is one of the 
clean decay modes of the KK graviton and in the region of interest, this could give 
an enhancement of the tail of the distribution. In Fig. [5] (top left panel) we have 
presented the invariant mass distributions of the photon pair for different choices of 
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Figure 4: Transverse momentum distribution of 73 [i.e, with minimum Pt\ (left 
panel) at the LHC for Ms = 3.5 TeV and d = 3. Rapidity distribution of 73 (right 
panel) at the LHC for Ms = 3.5 TeV and d = 3 m the region where P}^ G (750, 1250) 
GeV. 

Ms = (3.5, 4, 4.5) TeV fixing d = 3, while in the top right panel the same distribution 
is plotted for different choices of ci = 3, 4, 6 but a fixed Ms = 3.5 TeV. We find that 
the KK modes dominate over the SM contribution for larger values of invariant 
masses (say above 400 GeV for a given set of Ms and d values) of photon pairs 
leading to a significant enhancement of the signal over the background. We plot 
the factorisation scale dependence of invariant mass distributions of photon pairs in 
Fig. [6] (left panel) for different choices of fip, i-e., fip = (0.2, 2)Q. 

3. '-fZZ production: Invariant mass of Z boson pair is again a useful observable. 
We have done a similar analysis as we did for 77^ and use the same choice of 
factorisation scale and ADD model parameters. The invariant mass distributions 
are shown in the lower panels of Fig. [5] for different choices of Ms and d. We find 
that the invariant mass distributions of photon pairs in 77Z production and Z boson 
pairs here have similar qualitative behavior. In order to investigate the uncertainty 
resulting from the factorisation scale fip, in Fig. [6] (right panel), we have plotted the 
invariant mass distributions of Z pairs for different choices of ^f, ^-c, = (0.2, 2)Q. 

4. ZZZ production: We have classified the triple Z-bosons in such a way that 
P^^ > ^ > and for the P^' distribution we make the choice of factorisation 
scale /i_F = Pt'- Fig- El we have presented the transverse momentum distributions 
of Zi (left panel) and Z3 (middle panel) and rapidity distribution of Zi (right panel) 
for SM and SM+ADD with Ms = 3.5 TeV and d = 3. For the rapidity distribution, 
we have constrained 900 < P^^ < 1400 GeV. As in the case of 777, the P^^ 
distribution is similar to that of P^^ distribution. We have also shown the sensitivity 
of rapidity distribution to the factorisation scale hf by varying between hf = 0.2P^^ 
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Figure 5: Invariant mass distribution of the photon pair (top left panel) and Z 
boson pair (bottom left panel) at the LHC for d = 3 with different values of Ms- 
For Ms = 3.5 TeV with different values of d, the invariant mass distribution of the 
photon pair (right top panel) and Z pair (right bottom panel). 

and /i_F = 2P^\ In the central rapidity region we estimate the variation of the 
rapidity distribution with the factorisation scale and find that for fip = 0.2P^^ and 
fip = 2Py\ the variation is about 27.5 % and 8.9 % respectively with respect to 
fip = Pp^ ■ The rapidity distribution for Z2 is similar to that of Zi while Z3 is 
different. 

5. Total cross section: The total cross sections for various processes involving 
neutral triple gauge boson final states as a function of Ms for a fixed value oi d = 2 
are given in Fig. |8l The SM contributions that do not depend on ADD model 
parameter Ms appear as horizontal lines. 

So far in our numerical analysis we have put the UV cutoff A = M5 which is the 
conventional choice to do the phenomenology. The sensitivity of the choice of UV 
cutoff is presented in Fig.lHlfor P^^ distribution of 777 final state and also the invari- 
ant mass distribution of 77 pair of 77Z process by varying A = (0.9, 0.95, l)Ms. 
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Figure 6: Dependence of invariant mass distribution of the photon pair (left panel) 
on factorisation scale at the LHC for d = 3 and Ms = 3.5 TeV. Dependence of 
invariant mass distribution of Z boson pair (right panel) on factorisation scale at 
the LHC for d = 3 and Ms = 3.5 TeV. 

The cross section at P^^ = 1200 GeV varies between 10 - 24 % as we vary A = 
(0.9, 0.95)M5 as compared to A = Ms for the 777 process. Similarly the cross 
section for 77^' process aX Q = 2000 GeV the variation is 7 - 15 % in the same rage 
of A. 

The SM integrated K-factor for the 777 process is about 2.6 and for 77^ it 
is 1.5 and these are larger than uncertainties resulting from the factorization scale 
variations at LO [I6]. This warrants a full next to leading order QCD analysis which 
is beyond the scope of this paper. Such an analysis is reserved for future publication 

m- 

4 Conclusion 

In this paper, we have studied the neutral triple gauge boson production at the 
LHC in theories with large extra dimensions which are produced via the exchange 
of a tower of KK graviton taking into account the SM contributions. All the final 
state photons and Z-bosons are taken to be real. We have performed various checks 
on our analytical results, and the numerical predictions are obtained using a Monte 
Carlo code which allows us to implement various experimental cuts. Our code has 
also been tested to reproduce SM transverse momentum distribution of the largest 
transverse momentum for the 777 process. For the case in which the gauge bosons in 
the final state are identical we have presented the transverse momentum distribution 
by ordering the transverse momentum as > > Pj^. We find that P^ and P^ 
distributions are similar but the one for is different. The rapidity distributions 
are also presented. For the case where one of the gauge bosons in the final state is 
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Figure 7: Transverse momentum distribution of boson [i.e, with maximum Pq^ 
(left) and Z3 boson [i.e, with least Pt\ (middle) at the LHC for Ms = 3.5 TeV and 
d = 3. Rapidity distribution of that Zi boson (right) at the LHC for Ms = 3.5 TeV 
and = 3 in the region where P^' G (900, 1400) GeV. 



different we choose to use the invariant mass distribution of the identical di-bosons as 
it would be a better discriminator in the region of interest. We have also studied the 
dependency of the ADD model parameter M5 and the number of extra dimensions 
d keeping the UV scale A = M5. In addition we have reported the sensitivity of 
the choice of A by varying it from A = O.9M5 to O.95M5. We have also studied the 
dependence of our LO predictions on the factorisation scale. 
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Total Cross-section 




Figure 8: The total cross section for all the triple neutral gauge boson as a function of 
the extra dimension scale Ms for number of extra dimensions d = 2. The horizontal 
lines corresponds to the various SM contribution. 
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